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ABSTRACT 


When molecules of a low density diatomic gas strike a solid 
surface, both the translational energy and the internal energy modes 
of rotation and vibration will contribute to the energy exchange that 
occurs. 

The experimental apparatus described here uses the electron beam 
fluorescence detector to measure simultaneously both the rotational 
and the translational energy accommodation coefficient of room 
temperature nitrogen reflecting from a heated solid surface. Nickel, 
silver, gold and stainless steel surfaces were tested in a bakeable 
ultrahigh vacuum system built to provide a clean vacuum environment 
for control of the solid surface composition. 

Theoretical studies indicate that the accommodation coefficient 
of rotational energy should be less than that for translational energy, 
and this is confirmed by experimental results. 

Results show that for the pure metals - nickel, silver and gold, 
the translational accommodation coefficients a, are 0.5, 0.4 and 
0.2 respectively at 400°K; and the variation of Or with the increase 
in surface temperature, Te: is given approximately by the equation: 


4 


da,/dT, = -13.7278 x 10™" + 2.6042 x 107 /u 


where u is the ratio of the masses of the gas molecule and the solid 


atom. For the alloy stainless steel, the translational coefficient 
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decreases sharply with increase in surface temperature. 
Based on the above equation, an empirical relationship has been 


developed for translational accommodation coefficient a, in the form 


da, 
dr = qT, M + constant 


T 


where the constant of integration has been adjusted by fitting the 
equation to the best available experimental data. The resulting 
closed form solution for a, not only gives a good description of the 
results obtained in the present work but also predicts good results 
for a wide range of gas-surface combinations and conditions. 

The results also show that the rotational coefficient, Gp » is 
approximately 0.5 for stainless steel and 0.15 for nickel and gold. 
There is very little variation with surface temperature for these 
materials. For silver, the rotational coefficient increases from 
about 0.03 at 400°K to 0.2 at 700°K. 

In addition to being the only known direct method of measurement 
of the rotational acconmodation coefficient, the system used here 
offers an advantage over some previous methods of translational 
accommodation measurement in that there are few restrictions on 


solid surface temperature or composition. 
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T spectrometer rotational temperature measured by using 
the spectrometer 
Leate true rotational temperature 
Ly translation accommodation coefficient 
Op rotional accommodation coefficient 
a. internal accommodation coefficient 
) function of rotational quantum number and temperature 


u ratio of molecular weight of gas to the atomic weight of solid 


x14 


xosanortoaqe ott yd bewstv meod nowdosls to Mbtw 
nsowded sottixo sit woled mssd novtosis to eonst]eth 
4 bas A evadneds _ 
2owisiegneT 
26p insbront to swihsvsqmes nottstensyd 
asp bedaattey to siwfaveqmed nottslensts 
asp inabtont to swistsqmet Manofissor 
7 
asp betositer to suisveqnet fenoftsz01 a . 
swutsiaqnat Isnottstoy qt : 
ou ersqmsd sant we 
gftrau yd bevwesen swiIsvagmss lenortstoy “od amoit9 me 
rotanoritoege ott 
swtsreqmed Tenottsioy suit = 
oa sn" 
7 
dnataitseo sethdinmncis notanfen 7 


wave number 

wave number of rotational line (K',K'') = (3,2) in the 
(0,0) band of the first negative system 

density of gas 

solid angle subtended by the orifice at a point P 
constant defined in text 

Debye temperature of material of the solid 


angle noted on Figure 5 


ratio of specific heats 
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CHAPTER I 
INTRODUCTION 


At high temperatures and low densities, a gas in canonically 
distributed equilibrium follows the equipartition theorem. A diatomic 
gas in such an equilibrium has its internal energy divided between 
translational, rotational and vibrational degrees of motion in 
proportions which are kept in equilibrium by constant intermolecular 
collisions. This equilibrium can be disturbed in many ways. The 
passage of a shock wave, for example, would cause such a disturbance. 
From studies made on the passage of a shock wave through low density 
gases, it has been established that rotational and vibrational modes 
require many more collisions to come to equilibrium than the trans- 
lational mode does. 

Since energy exchange between a solid surface and a gas is also 
accomplished by collisions, it can be expected that the energy transfer 
to the internal modes is less than that to the translational mode.. The 
accommodation coefficient provides a measure of the energy exchanged to 
each mode. 

As the energy exchange between a solid surface and a gas obviously 
depends on the mechanism of surface atom-gas molecule interaction, many 
attempts have been made to predict theoretically, the results of such 
collisions. Thus, Maxwell, in extending to low pressure Poiseuille's law for 
capillary flow, suggested that a fraction f of the molecules colliding 


with a surface would be adsorbed and be diffusely reflected, while 
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(1-f) of the molecules would reflect specularly. Only the momentum 
history of the adsorbed molecules was considered. Nothing was said 
about their energies. M. Smoluchowski (1) used this idea and 
introduced a similar fraction f, which were adsorbed and reflected 
at the temperature of the surface. He showed this fraction f to be 
the same as the acconmodation coefficient defined by Knudsen. 
Stickney's (2) experiments measuring normal momentum transfer of a 
molecular beam directed against a metal surface, led him to the 
conclusion that a fraction are fully accommodated and reflect 
diffusely with a velocity distribution corresponding to the surface 
temperature, while the remainder, experience no energy accommodation 
but are diffusely reflected. No special precautions were taken to 
remove adsorbed gas layers from the solid surface in these experiments 
and diffuse reflection can always be expected from gas covered surfaces. 
A review of Knudsen's work can be found in Loeb's book (3) on 
the Kinetic Theory of Gases. The following account is derived from 
this source. Knudsen employed the kinetic theory of gases rather 
adroitly in his investigation of heat transfer at very low densities 
between two parallel plates separated by a distance less than the 
mean free path of the gas. He made the assumption that both the 
plates were "molecularly rough", meaning that gas molecules which 


collided with the plates left them at the temperature of the plates. 
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From kinetic theory, he obtained an expression for the energy 


transferred to the cooler plate as 
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where k ayy. y is the ratio of the specific heats cp/cy ; 


E 


T is the absolute temperature of the gas, and Po is the density in 
grams per cubic centimeter of the gas at 273°K and a pressure of one 
dyne per square centimeter. Ty and Ts, are the temperatures of the 
plates and P is the pressure of the gas. Since E itself is difficult 
to measure because of radiation and residual gas conduction losses, 


measuring E at two pressures for a given Ty - To gives the ratio 
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aT Mi 1.3 
where M is the molecular weight of the gas. 

Experimental results gave lower values for this quantity compared 
to those obtained from equation 1.3. At this juncture, Knudsen 
questioned the assumption of molecular roughness. He decided that 
the energy transfer between the gas and the surface was not complete 
as had been assumed. This led to the definition of accommodation 
coefficients. Knudsen defined the coefficient based on temperatures, 
which assumes Maxwellian equilibrium distribution at rest for the 
molecules. Much later, he conducted an experiment in which he 


obtained the translational acconmodation by radiometric forces and 
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the total accommodation by the conduction of heat. In developing 
the theory for this experiment he assumed that the rotational energy 
of the molecules corresponded to that of a canonical energy 
distribution at the temperature of the solid surfaces. The results 
obtained for hydrogen and helium on a platinum strip, bright on one 
side and roughened on the other showed that the accommodation 
coefficients for the translational mode and for the internal mode 
were equal. 

In modern literature, the coefficient of acconmodation is 


defined as 
. ~ EE.) 1.4 


where E. is the energy of the incident gas, E. is the energy of the 
reflected gas, and E. is the ideal energy if the gas is reflected 
at the temperature of the surface. This definition reduces to 
Knudsen's original definition of a based on temperatures, if both 
the incident and the reflected gas molecules have a Maxwellian velocity 
distribution. 

In the present work, two accommodation coefficients have been 
used: the translational accommodation coefficient Or » and the 
rotational accommodation coefficient op. These are defined as 


follows: 
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Logan and Stickney's (5) hardcube theory uses a simplistic 
model, proposed mainly to reduce the mathematical complexities of the 
more realistic models such as those of Goodman (4), and to take into 
consideration the thermal motion of the surface atoms. These are 
assumed to be hardcubes which have vibrational motion perpendicular to 
the surface. The interaction potential between these cubes and a 
gas molecule is a step function repulsive force which does not allow 
the gas molecule to spend any time on the surface. No trapping of the 
gas molecules is possible. 

The results from such a model for molecular beams show one 
feature of interest: as the surface temperature increases the 
intensity maximum of the reflected molecule shifts towards the surface 
perpendicular, in agreement with experimental results cited by the 
authors. 

Logan and Keck's (6) softcube theory is a logical extension of 
Logan and Stickney's (5) theory and assumes the potential between the 
solid cube oscillator and the gas molecule to be exponential. The 
theory was proposed mainly to provide better agreement with experimental 
results than the hardcube theory does. It also takes into consideration 
the trapping of gas molecules. The agreement given by the hardcube 
theory for the behaviour of the intensity maximum with temperature 


comes out of this theory also. 
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Among theoretical papers dealing with the internal mode of 
energies, those of Oman (7), Feuer (8) and Marsh (9) are of interest. 
Oman has performed a computer simulation of the collisions 

between an incident gas of homonuclear molecules and a solid surface 

made up of atoms in a lattice. Again, the lattice has been taken 

to be at O°K. The interaction potential is 6 - 12 Lennard-Jones 

type. The results show that the rotational energy exchange varies 

linearly with the translational energy of the incident molecules. 
Feuer (8) has made a theoretical study based on-a quantum 

mechanics model assuming a single phonon transition in the solid. 

The results indicate that the rotational accommodation coefficient 

is less than the translational accommodation coefficient. Ina later 

paper, Feuer and Osburn (10) have calculated the two accommodation 

coefficients of para and ortho hydrogen on platinum at two temperatures, 

one above and one below the Debye temperature of the solid. Table 


1. below shows their results and is reproduced from (10). 


Table 1 


Theoretical values of translational and internal accommodation coefficients 


a para a. para/a ortho 


338 0.337 0.0344 
150 Wasi 0.0423 
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The internal energy coefficient is very much less than the translational 
coefficient for both para and ortho hydrogen. The results given in 
Table 1. were obtained by considering only single phonon transition 
of the solid atoms. Two phonon transitions have been dealt with by 
Allen and Feuer (11) and they show that it contributes very little 
towards accommodation. For helium on platinum, for example, 
accommodation due to two phonon transitions is 1.33% of that due to 
single phonon transitions. Therefore the results considering only 
Single phonon transitions should be good. 

Marsh (9) has devised a gas-surface interaction model based 
on physical adsorption. Translational and internal energy 
accommodation coefficients are obtained from equations similar to 
gas-phase relaxation equations. The resulting expressions contain 
the relevant heats of adsorption. Results for diatomic gases, 
pertinent to this work, have not been obtained due to paucity of 
information about the heats of adsorption. 

Experiments to determine thermal accommodation coefficient can 
be classified into two main categories: the first, in which the total 
thermal accommodation is obtained, and the second, in which both 
translational and internal coefficients are obtained. 

Among the many methods belonging to the first category, Knudsen's 
low pressure method and the temperature jump method are perhaps the 


most popular. These methods are well known and have been documented 
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by various authors, including Wachman (12). Apart from the fact that 
these methods give only the total thermal accommodation, they are also 
limited in the temperature and composition of the solid surface that 
can be used. 

Translational accommodation coefficient can be determined from 
molecular beam experiments. |The method usually consists of shooting 
a monoenergetic beam of molecules at a given angle to the test surface, 
and then determining the number density and velocity of the reflected 
molecules at various angles with respect to the surface. Normal 
momentum accommodation is obtained from this. If the reflected 
molecules are assumed to have a Maxwellian velocity distribution, 
then the translational accommodation coefficient can also be determined. 
Accommodation of other modes is not determined by this method. 

Of more interest and concern to the present work are the 
experiments belonging to the second category. Since there were no 
direct methods for measuring the internal energies until the electron 
beam fluorescence technique was developed, most of the experiments in 
the second category measured the internal energy accommodation by 
subtracting the translational accommodation from the total thermal 
accommodation. 

Sasaki, Taku and Mitani (13) used a torsion balance to determine 
a The test surface, in the form of a strip which can be electrically 


r 
heated, is attached to one wing of the balance. A horizontal molecular 
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beam is shot at this strip. The translational coefficient is a 
function of the angles to which the balance rotates, depending on 

the temperature of the strip. The total thermal coefficient is 
obtained by measuring the heat loss from the strip when the strip is 
at the higher temperature (Knudsen's low pressure method). They 
obtained the values for dr = 0.393 and Ap = 0.314 for a nickel strip 
at 518°C with nitrogen at 31.3°C. 

In a later paper, Sasaki and Mitani (14) used a circular cross 
sectioned filament, shaped like an inverted V. They found the 
translational accommodation by relating it to the period of half 
decrement when the filament is set vibrating, with the translational 
motion of the molecules providing a damping mechanism to impede 
vibration of the filament. Again, the total thermal coefficient was 
found by Knudsen's low pressure method. They obtained ay = 0.48 
and Cp = 0.10 for tungsten at 300°C with hydrogen at 23°C. 

Variations and modifications of the torsion balance method have 
been used by Stickney (2), Carpenter, Humphries and Mair (15) and many 
others, but mostly to obtain normal momentum coefficient or translational 
energy accommodation coefficient. 

Schafer and Riggert (16) used a method claimed to be more accurate 
than that used by Eucken and Krome (17) to obtain both translational 
and internal energy accommodation coefficients. They heated a band 


of the test material in a Knudsen cell. A filament of the same 
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material fixed parallel to the band picked up heat from the reflected 
molecules. From this they were able to measure the thermal 
accommodation of the filament. They also determined the thermal 
accommodation coefficient of the filament alone by the low pressure 
Knudsen method. 

Now, if all the modes of energy were equally accommodated, 
thermal accommodation found for the filament independently would be 
the same as that found when it is heated by the reflected molecules 
from the band. On the other hand, systematic differences in the two 
values can be related to partial accommodation coefficients. 

Their results show that for nitrogen on gold, ou decreases 
from .89 to .86 as de increases from 293°K to 389°K, and Op decreases 
Frohne 7B to .54 for the same temperature range. 

A technique for direct measurement of rotational energy 
distribution of nitrogen molecules after reflection from a target surface 
was developed by Marsden (18), using electron beam excited fluorescence 

of light. Fluorescence from the beam at a point in front of the 
target was made up of contribution from molecules reflected directly 
from the target surface and from molecules of the background gas. 

The contribution from the background gas could be determined separately 
by retracting the target and this was then subtracted from the first 


measurement to obtain the contribution due to reflected molecules alone. 
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Results for a silver target showed that translational and rotational 
accommodation coefficients were equal and had a value of 0.85. _—siIt 
may be noted that the rotational energy distribution of the reflected 
molecules was found to correspond to that of a gas in thermal equilibrium. 
These results, were for a gas at higher temperatures than the surface. 

The apparatus described in this thesis uses the electron beam 
fluorescence detector to provide simultaneous direct measurement of 
translational and rotational accommodation coefficients for nitrogen 
on a solid surface. It is an extension of the technique used by 
Marsden. 

Two parameters were deemed to be of importance: the temperature 
of the surface Ler and the ratio u of the molecular weight of the gas 
to the atomic weight of the surface material. These parameters were 
varied during the experiment, keeping the incident gas nitrogen at 
room temperature. Measurements were made for a temperature range of 


de = 300°K to 900°K, and a u range of 0.142 to 0.483. 
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CHAPTER II 
EXPERIMENTAL APPARATUS 


2.1 Vacuum System 

The stainless steel vacuum System was constructed using inert 
gas welding technique for permanent joints and Varian Conflat type 
flanges with bakeable copper gaskets for demountable joints. An oven 
made of commercially available insulating material and covered with 
aluminum foil, enclosed the system and generated bakeout temperatures 
of about 300°C. A schematic of the system is shown in Figure 1. The 
different chambers have been labelled A, B and C to facilitate 
explanation. The top chamber A, is connected to chamber B by a 2.08nm. 
diameter orifice, as well as by a 37mm.diameter valved bypass. The 
bypass helps in reducing pressure in chamber A during bakeout. Chamber 
B sits atop chamber C and is separated from it by a large passage which 
presents virtually no obstruction to the flow of gas. The pressure 


in the two chambers are equal. 


2.2 Chamber A 

The gas-surface interaction takes place in chamber A. The gas 
introduction system on the left side of Figure | consists of a big 
plenum chamber from which gas is led into chamber A via a cold trap 
and two leak valves. A commercial gas bottle and a pump-stand are 


connected to the plenum chamber with the help of a tee-joint and 
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appropriately located stop-cocks. The plenum chamber contains gas 
at high pressure and provides constant pressure upstream of the leak 
valves. The cold trap eliminates any moisture present in the gas and 
the leak valves are used to meter the flow of gas into the system. 
Chamber A is watercooled to keep the temperature of the gas constant. 
The test surface is made of 0.005 or 0.01 inch thick, 7mm. wide 
Sheet metal. Figure 2 shows the details. The target is U shaped 
lengthwise, the bottom flattened out to a plane surface about 10mm. 
long. The arms of the U are attached to the pinched ends of two 
heavy stainless steel tubes by screws and nuts. Both tubes are 
attached to the flange on top of chamber A, one of them being insulated 
from it. The location and other details and dimensions of the target 


are shown in Figure 2. 


2.3 Chamber B 

The roughing pump system is attached to chamber B on the right 
side in Figure 1. The pump system consists of a 4 litres per second 
rotary pump and a 100 litres per second oil diffusion pump. It is 
attached to chamber B by 37mm. stainless steel tubing. A cold trap 
is present in the connecting circuit. Rough pumping takes the vacuum 
system to about 1 x 107° torr pressure. 

Chamber B also contains the electronbeam. A Phillips 902-564 


television type electron gun to the left of chamber B produces a well 
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Figure 2. Target Mounting Details. 
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collimated, electrostatically focussed 1.5mm. diameter beam. The 
mechanism by which the gun is attached to chamber B facilitates 
adjustment of the location and direction of the beam. As shown in 
Figure 3, the beam passes through four 2.5mm. diameter holes, which 
locate it directly below the 2.08m. orifice between chambers A and B. 
The power for the gun is supplied by an N.J.E. 0-30KV, 0-2.5ma. power 
supply. At normal operating filament current of 2.0ma and accelerating 
voltage of 10KV, about 1.1lma.beam current is obtained. This is picked 
up by the receiver on the right side of chamber B and is monitored by 
a digital milliammeter. 

Some stratigically placed light baffles in chamber B prevent 
stray light from reaching the light sensing devices located outside 


a viewing port on chamber B. 


2.4 Chamber C 

Chamber C is the pumping chamber. It is large enough to prevent 
sudden fluctuations in pressure. A 100 litre per second triode ion 
pump and a titanium sublimation pump provide the pumping. To prevent 
heating of the gas by the sublimation pump , there is a liquid nitrogen 
cooled inner shell in chamber C. This shell also provides a cool 
surface upon which the sublimating titanium molecules can settle. Only 
the ion pump was used while baking and for the final pump down. It 
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provided a pressure of 7 x 10 ~ torr. The sublimation pump was not 
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used during any part of the experiment as it proved incapable of 


providing constant pumping speed. 


2.5 Pressure Measurement 

A nude ultra high vacuum Varian ionization gauge attached to 
chamber C was used for the measurement of background pressure. As 
will be evident from the theory of this experiment, the background 
pressure needs only to be maintained constant. Since measurement 
of absolute pressure was not necessary, no accurate calibration of 
the ionization gauge was needed. The pumping current of the ion 


pump also provided a rough check on the total pressure. 


2.6 Temperature Measurement 

The target was heated ohmically by a 50 ampere 40 volt direct 
current power supply, connected to the copper feed through at the top 
of chamber A. A Chromel-Alumel thermocouple was spot welded to the 
target to measure its temperature. Calibration of the thermocouple 
was carried out outside the system and results were within 71°C of the 
values given in pertinent tables. The temperature gradient along the 
length of the target was also found by scanning with a radiation 
thermometer and it was ascertained that no cold spot was formed at the 
attachment point of the thermocouple. 

The temperature of the test gas introduced into chamber A was 


assumed to be the same as that of the wall of chamber A. Two Chromel - 
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Alumel thermocouples attached to the wall measured the temperature. 
Any tendency for the wall to heat up due to radiated heat from the 
target and the leads was compensated for by adjusting the flow of 
cooling water. The walls were maintained at room temperature at all 


times. 


Zip SOpties 

The schematic of the optical system is shown in Figure 3. A 
cylindrical concave mirror, with its center of curvature coincident with 
the electron beam, reflects the fluorescent light to a 150mm. focal 
length double convex lens. The lens, which is outside the vacuum 
system in front of a viewing port in chamber B, focusses the image of 
the beam onto the entrance slit of the f/6 Czerny-Turner grating 
Spectrometer. The spectrometer was placed on its side so that the 
slits were parallel to the beam. The intensity of the resolved 
Spectrum lines was measured by an EMI 9502S photomultiplier. 

The power for the photomultiplier was supplied by a Fluke model 
410B 0-2.1 KV supply. It requlates output voltage within 0.001%. At 
normal operating conditions of -1.5 KV, the photomultiplier has a dark 
current of 4 x 1071? amperes. 

The output from the photomultiplier was measured by a Keithly 
610 B Electrometer and recorded by a Honeywell Electronik 19 chart 


recorder. 
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CHAPTER III 
EXPERIMENTAL PROCEDURE 


3.1 Surface Target Preparation 

The target was cut out to the required size from 0.005 or 0.010 
inch thick high purity metal foil. One face was cleaned with acetone 
and distilled water, and thermocouple wires were spot welded on it 
with a dial-a-weld telephone type welder. The target was then bent 
to shape and attached to the electrical leads as shown in Figure 2. 
Both the target and the leads were thoroughly cleaned with acetone 
and distilled water and dried with hot air. The whole assembly was 
then mounted onto chamber A, taking care not to leave greasy spots or 
fingerprints. |The vacuum system was roughed, baked and pumped down 
to 7 x 107? torr pressure. 

The target was heated to a high temperature (900°K in the case of 
nickel and gold, 750°K in the case of silver) at this pressure. 


° torr due to outgassing, 


Initially, the pressure went up to about 6 x 10° 
but after an interval of about 2 hours, it settled down again to 
Ue ee torr, in which condition it was maintained for a period of 
twelve hours. 


The stainless steel specimen was not cleaned in as thorough a 


fashion, since experiments on it were carried out while assessing the 


experimental method. 
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3.2 Introduction of Test Gas 

The plenum chamber, the cold trap, and the manual leak valve were 
first pumped down to about 1 x 1076 torr pressure. Nitrogen from the 
bottle was then introduced into them at about 5 psig. They were 
again pumped down. This process of flushing the system was carried 
out at least ten times to insure that the introduction system contained 
only nitrogen. 

The plenum chamber was again filled with nitrogen to a pressure 
of about 5 psig. Every effort was made to keep this pressure constant 
by operating the bottle regulator whenever the gauges. on the regulator 
showed a decrease. At this high upstream pressure of 5 psig, the servo 
driven needle valve failed to perform satisfactorily. It was therefore 
left open permanently. The manual needle leak valve was used to 
introduce gas into the main vacuum system and to control the background 
pressure. 

Before any of the experiments were begun, a partial pressure gauge 
was attached to chamber C and the partial pressure of nitrogen at various 
total pressures was determined. Figure 4 is a graph of the results 
obtained. Tests were repeated a number of times and it was ascertained 
that the results shown on the graph did not vary. 

After the target had been kept hot for 12 hours, the test gas was 
introduced into the system, and with the aid of the manual leak valve, 
the pressure as read by the ionization guage in chamber C was steadied 


to a constant value between 2 x 1076 torr to 4 x 107° CONT: 
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3.3 Determination of Accommodation Coefficients 

After the target had been cleaned by heating, maintaining it at 
the same high temperature, the R branch of the N, (0,0) band emission 
was slowly scanned with the spectrometer. The results were recorded. 
Scans were carried out with the bypass valve between chambers A and 
B both open and closed. Measurements were made at a series of target 
temperatures varying from bake out temperature to room temperature in 
100°K intervals. 

The speed of scanning, the pressure measured by the ionization guage, 
the beam current, and the temperature of the incident gas were all kept 
constant during the entire duration of the experiment. 

The experiment was repeated a number of times to test for 


repeatability. 
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CHAPTER IV 
THEORY 


4.1 Geometrical Considerations 

A heated solid surface made up of the material to be tested is 
located in chamber A opposite a small orifice leading to chamber B as 
shown in Figure 1. 


The pressure in chamber A is of the order of ile 


torr, and the 
target surface is placed in such a way that substantially all the 
molecules reaching the electron beam from chamber A have come directly 
from the target surface with no intermediate collisions. The pressure 
in chamber B is kept as low as possible, typically 10° torr. 

The valved by-pass connecting chambers A and B allows for pumping 
to outgas impurities that may be absorbed on the system walls. This 
by-pass is also used to reduce the pressure in A and thereby reduce the 
contribution of molecules coming through the orifice to a small fraction 
of the background density so that the contribution of the background 
gas in chamber B can be determined. 

The number density of gas at a point P is n = n, + Nb where n. 
is the contribution to number density by the molecules reflected from 
the target and issuing through the orifice, and ny is the contribution of 
the background molecules in chamber B. 

Since the gas in chamber A is in equilibrium and at rest, the 


velocity distribution is Maxwellian. It will be assumed that the 


reflected molecules come from a gas in equilibrium and at rest with a 
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temperature she and number density n “ above the orifice. Consider a 
small area dA sce point P subtending an angle dO at centre of orifice, 
at an angle 6 to the vertical, as shown in Figure 5. 

The number of reflected molecules in solid angle dw passing 
through dA with a velocity c in time dt is: 


N(c) = Ge cos@ dcdt dw dA/cosé pe 


Where dw is the solid angle subtended by the orifice at point P and 

fe is the aeons distribution function of the reflected molecules. 
Now dw = Aucos 8 and dA = ae wd6@, where w is the width of the beam 
which is serie through the slit of the spectometer, and A' is the 
orifice area. 


The number passing through dA in time dt with velocity c is, 


N(c) = cF at cos*ede dedt 4.2 


The number density of reflected molecules in dA is the number 


passing through dA divided by cdAdt, 


el an A; cos"6 dc hae) 
y 


The number density due to reflected molecules of all velocities 
is obtained by integrating equation (4.3), 


i ot, a cos26 dc = Bet ens 4.4 
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The number density in dA due to background molecules is n Zz: 


g | 
The average number density in the length of beam viewed by the 


optical system is obtained by taking the integral 


6, Nn 
Oar Bra | 3 dw 
i, IZ > cos 8 + np (1 - 7) | A 


where + 8) marks the ends of the portion of electron beam viewed. 


' : 6 
4 A Sanya G 
(n.-ng) ne joy + GS | 
n= y + np 
tan 1/4 + 8 /2 
'" tan a/4 ~ 8,72 4.5 
This can be written as, 
n=. neat (1 - 2) Np 4.6 


The constant & depends only on system geometry. 

When the bypass between A and B is closed, the orifice is the only 
connection between the two chambers. Let its conductance be C,. Also 
let the speed of the triode ion pump be S, and the conductance of the 
bypass be C,. eee, chamber A is related to that in chamber B 
by the relation Pa = ( Pp when bypass is closed, and since the 


S+C 
temperature in the two chambers is same, Ny = (——)n,. From conserva- 
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substituting in (4.6) 


S+C 
7 1 
n = Q( e ) eae np + (1-2)n, 4.7 
and similarly when the bypass is open, 
S+C,+C, 
i = rome oe) OB bu (1-2)n, 4.8 


4.2 Electron Beam Measurements 


The high energy electron beam excites and ionizes the neutral 


] 


nitrogen molecules, which are initially in N, X E ground state, to 


Z 
an upper state he The light emitted comes almost entirely 


(19) from the first negative band due to the transition No" Bor" = 


+ (2 + 


Ny X oF of the nitrogen ion. If the doublets are not resolved, the 


transition becomes 'Z > 'Z and the applicable selection rule AK = +1, 
gives rise to two branches of emission P and R. The P branch has a 
vertex towards red on the Fortrat diagram and forms a bandhead. The 
R branch has peaks occurring at ever widening intervals towards the 
ultra-violet. A scan of the R branch of the N5(0,0) band is shown in 
Figure 6. The two to one alternation in intensity between odd and 
even numbered lines is a characteristic of nitrogen. The excitation 
emission mechanism is given in more detail in the Appendix. 

The number density is determined by a summation of all the light 
output from one branch of the N5(0,0) band or alternatively by measuring 


the total light output from the unresolved band. Light intensity is 


given by the equation, I = nkp, where n is number density, k is constant 
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depending on the optical system and beam current and is a factor 
depending on the spectrum lines being observed. If the rotational 
spectrum is resolved o(K' Tp) depends on rotational quantum number K' 

and rotational temperature Tp» n and k being constant. If the spectrum 
is unresolved and all light from the band is observed I = nk®, where 

6 = i o(K',Tp). ® is a constant, and I is then proportional to n. 

; rte distribution of rotational energy, and the rotational temperature 

of the gas, can be obtained directly from a measurement of rotational 


line intensities (19). Rotational energy in a diatomic gas has the 


distribution (32) given by equation (4.9), 


Ny = aT (2341 )e7PM© 3( 31) /kT 4.9 
where Ny = number of molecules in rotational energy level J 
N = total number of molecules 
B = the rotational constant for the molecule 
h = Planck's constant 
k = Boltzmann's constant 
c = speed of light 


In the case of the N5(0,0) band of the nitrogen spectrum the line 
intensity distribution is given (19) by equation (4.10), 
- thal Oa 
Teese sae rete Kage eet Mxql6l “ze (K +1)hc/kTp 4.10 


0 
where K is the rotational quantum number apart from spin (J = K+l1/2). 
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A single prime refers to the upper electronic state and double prime 
refers to the lower electronic state. 
i] 1 he t hc 1 
(K'+1)exp[-2 B (K'+1);=—] + K' exp[2B — K'] 
kT, kT, 


NS 


(2K' + 1) fai 


Ea Bs 


[G] is a function of Tp» 
Xq is a constant related to overall intensity, and 


V/V is approximately constant for a given band. 


Equation (4.10) is derived from equation (4.9) by going through the 
process of excitation (and ionization) of the nitrogen molecules by the 
electron beam, followed by a spontaneous decay to the ground state of 
the ion with emission of light. The temperature Tp and rotational 
constant B appearing in both equations are the same. 

A plot of Log[I(K')/(K'+K''+1) [@]xq(v/vg)"] against K'(K'+1) will 
result in a straight line with a slope ~Bhc/kT, if the gas has a equi- 
librium rotational energy distribution. 

Thus the rotational temperature Tp can be determined from a 
measurement of spectrum line intensities. 

Figure 7 shows some plots of measured line intensities used to 


determine rotational temperature. 


4.3 Translational Accommodation Coefficients 
The translational energy acconmodation coefficient, Oy» has been 


defined in equation 1.5 and can be written as, 
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GAS INTRODUCED AT 295°K 
Tr =290°K (measured background) 
Ts =550°K 


Tr, =402°K (measured reflected) 
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Figure 7. Typical Line Intensities in the Background and 


Reflected Molecules. 
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On Shs Alan al hese Cheeni 1) 4.12 
il Lins T. S T; 


where Ty is the temperature of the incident gas, Ty is the temperature 
of the reflected gas and il is the temperature of the solid surface. 
The temperature ratio Ty /T+ is determined experimentally as follows: 


ahead | 
When the bypass is open, the number density measured at the electron 


beam is, 
S+C,+C, 
ny) = 2 ( 0 #C, ees Np ae eae 1 (1-2)n, 4.13 
and when the bypass is closed the number ee measured is, 
S+C, 
No = Q (= eg) Np v1, + (] -2)Nnp 4.14 


Measurements were made with the target temperature the same as that of 


the incident gas ie = Ty = Ty ), and at a higher temperature resulting 


r j 
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Since I = kn® and since k and @ are constants, 


17 
where I = = ae and subscripts 1 and 2 refer to the bypass open 
K's 
or closed respectively. 


From equation (4.12) 


Or = 1 if 
ae ip Dr, 
TIT; > 4.18 
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4.4 Rotational Accommodation Coefficient 

Line intensities of the R-branch of N,(0,0) band for the molecules 
reflected from a heated surface are obtained in this experiment by 
subtracting the line intensities scanned with the bypass valve between 
chambers A and B open from those when it is closed. The rotational 
temperature of the reflected molecules is determined from these measured 
line intensities. 

The number density n from equation 4.6 is 


n= On, + Cit = Q)np 


The intensity of a measured spectrum line corresponding to 


rotational level K' can be written as 


Lhe Qn. ko(K', TR!) + (| = 2)n, ko(K', We 4.19 


where, as before, k depends on the geometry of the optical system and 
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@ is a function of K' and the temperature of the gas. 


Equation 4.19 can be written as 


I(K') = QIn(K'sTp.) + (1-2) In (K",T 4.20 


R.! 
F 


where In(K' 51, ) is the contribution to line intensity at K' of the 
re 


R.) 


reflected gas molecules at a rotational temperature Tp » and I,(K'T 
r j 


is the contribution of the background molecules at temperature Tp 5 
i 


As indicated before, the R-branch is scanned with the bypass valve 
open and closed, keeping all the other variables constant. The gas 
number densities in the two cases are examined below. 

When the bypass valve is open, the number density n is given by 
equation 4.8. Substituting S = 100 liters per second, Cy = 20 liters 
per second, C, = 0.4 liters per second and 2 = 0.001876 (obtained by 
substituting appropriate values in equation 4.5), equation 4.8 becomes 


n = 0.011072 ng ey + 0.998124 ng 4.21 


Ty 


where the first part of the right hand side of the equation is the 

contribution of the reflected molecules at a translational temperature 

Ty , and the second part is the contribution of the background molecules. 
r 


In the present experiments, I was always greater than or equal to 
r 


T. because the surface target temperature was always greater than or 


Me 


equal to the temperature of the incident gas. In view of this, the 
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Equation 4.21 reduces to 


n = 0.011072 n, + 0.998124 n 4.22 


B B 
which reveals that the contribution of the reflected molecules to the 
total number density is about 1% while the background molecules contribute 
99%. 
At this juncture an approximation will be made that 
n= Np 4.23 
when the bypass valve is open. This approximation will be eliminated 
by a refining process, shown later. 
With this approximation equation 4.20 reduces to 
I(K') = In(K' T 4.24 


R.) 
F 


when the bypass valve is open. 

Since these line intensities are from background gas molecules known 
to be at 295°K, instead of taking I(K') directly from the spectrographic 
scan, the corresponding theoretical values are taken to eliminate 
unnecessary scatter in the results. 

The slope of the line used to determine the rotational temperature 


is Bhc/kT, , and for Tp = 295°K has a value of 0.004237 on the common 
i 


log scale. 

Using the negative of this slope, an equation for a straight line 
is written as Y = -0.004237 x+C 
where ge log U(K*) (4K! 41) [6] 29") and % = K'(K’ 41). 


C is an arbitrary constant, and for the present can be taken as 0. 
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From this equation x and y are calculated for various values of K'. The 
calculated values are given in Table 2. For even values of K', I(K') 


has half the value given by the equation as explained in the appendix. 


Table 2 


Theoretical values of line intensities at room 
temperature 


log (k'+k"'+1)[6] (2)! 


In column 4 of Table 2, [G] is calculated from equation 4.11 and values 
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of v are taken from Child's (20) paper. To distinguish calculated 
values from measured values of intensities, let the values given in 


Table 2 be called Lanne Ke Je 


Values of Le apletk') from K' = 3 to 15 are summed up. This sum 
is equal to 60.9535. Values of I(K') from the experimental scan of 


the R-branch with the bypass valve open are also added to obtain 
15 


ZX I(K'). In the first approximation, I,(K') is taken to be 
K'=3 
15 
I,(K") = aS CK) SETI 
K'=3 
Now, when the bypass valve is closed, the number density n is given 


table Kk’ )/60.9535 Anco 


by equation 4.7. 
Substituting values for S, Cy» and 2, 


n = 0.470n., VT, /T+ + 0.998124n, 4.26 
j Yr 


where, again, the first part of the right hand side is the contribution 


of the reflected gas molecules, and the other part is that of the back- 


ground gas molecules. In the present experiments, the maximum value 
that ils) reached was 537°K corresponding to Oe = 0.4 at bs = 900°K for 
r 

gold. This gives vT- /T- = 0.7425. Equation 4.26 becomes 


T; a 


n = 0.3496n, ot 0.998124n, 4.27 


It is of interest to note here that the contribution of the reflected 
molecules to the total number density varies between 0.3496n, and 


0.4708n, corresponding to Tr = 537°K and Tr = 295°K, 
ie r 
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Examining equations 4.26 and 4.20, values of the line intensities, 
I,(K'), due to the reflected molecules, can be obtained by subtracting 
In(K') (given by equation 4.25) from the line intensities obtained by 
scanning the R-branch with the valve closed. The result can be written 
as 

QI,(K') = 1(K') - Ip(K') 4.28 
A computer program was written to fit a least squares line to points 


given by plotting Vog@Ip(K')/(K'#K' 41) [612 )* against K'(K'+1). 
) 


This program is given in the appendix. From the negative of the 
Slope of this line, the rotational temperature of the reflected gas 


molecules Tp is obtained. 

As was none an approximation was made in writing equations 4.23 
and 4.24. The value of TR found above contains this approximation. 
A refining procedure to correct this follows. 


Substituting the value of the translational temperature i; in 
ie 
equation 4.21, the total light contribution of the reflected gas molecules 


can be written as 


15 
Ae = 1(K') 
srry Ot x 0.011072 vT-/T,_ 
- 
3 R 0.011072 72 + 0.998124 i, ft 
i 4.29 
and 
15 3 115 15 
zr I,(K") ese YUN Ag ge I,(K") 4,30 
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Values of ZI, (K') given by equation 4.30 are used to find the 
corresponding theoretical line intensities as was done in Table 2. 

These theoretical values are substituted into equation 4.28 to obtain 
I,(K") and subsequently the rotational temperature Tp from the computer 


\e 
program. 


In practice, however, it was found that the approximate equation 


4.24 for I,(K'); corrected as per Table 2, resulted in values of Tp 
. 


which were within 0.06% of the value found after refining. 
Once the rotational temperature Tp is found, the rotational 
r 

accommodation coefficient is given by equation 1.6 

elemtTpy = Tee eds (TAN TECD 1.6 

Yi 1 i 

The pumping speed S and the value of the conductance Cy are known only 
‘approximately, but this does not affect the accuracy very much. In 
equation 4.21, if the value of S is taken to be 50 liters per second 


instead of 100 liters per second 


t 


n= 0.00776n, V1 /T+ + 0.998124n, and 
Hl r 


Hi Tae is taken to be unity as was done to obtain equation 4.22, then 
n= 0.00776n, ip 0.998124n, 

which, when compared to equation 4.22, gives better justification to 

the assumption made in writing n = Np. It is doubtful that the value 


of S would exceed the 100 liters per second specified by the manufacturers. 
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Similarly, C, has been assumed to be 20 liters per second based on 
the geometrical dimensions of the 37 mm. inch tubing bypass connecting 
chambers A and B. Substituting two values which differ by ten liters 
per second from the value assumed for Co in equation 4.8, 

n = 0.0199In, + 0.998124n, 
for Co = 10 liters per second if Ty /Ty is unity, and 


1 ahr 
n = 0.00804/7n, m2 0.998124n, 


for Cy 30 liters per second. 

Comparison of these two equations to 4.22, which has Cy = 20 liters 
per second, indicates that the contribution of the reflected molecules 
(the first part of the right hand side) is altered by less than 1%. 


Thus, the values of S and Co» which are known only approximately, 


do not affect the accuracy of the method very much. 
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CHAPTER V 
RESULTS AND DISCUSSION 


5.1 Rotational Accommodation Coefficient 

Figure 8 shows the variation of rotational accommodation coefficient 
Op with target surface temperature is for nickel, silver and gold. The 
value of Op is about 0.15 for nickel and gold. There appears to be 
very little variation with T,. For silver, however, Ap varies from 
0.03 at 400°K to 0.2 at 700°K. 

Results of Sasaki, Taku and Mitani (13) show that at T, = 791°K and 
Use fy 304. 32K. Ap = 0.314 for nickel and nitrogen. Results obtained 
in the present work for the same temperature and surface give Op = O.19, 
which is about half in magnitude. The difference is probably due to 
the condition of the test surfaces. It is well known that accommodation 
coefficients are sensitive to surface condition. For example, it was 
noticed during the present experiments that stainless steel, which will 
have an oxide covered surface, exhibited high values of Gp (Figure 9). 
Sasaki, Taku and Mitani (13) have not given enough information in their 
paper to draw any conclusions about the condition or their solid surface 
or the reliability of their measurements. No attempt appears to have 
been made to clean the surface, other than working at high temperatures. 
Sasaki and Mitani (14) in a later paper give Gp = 0.10 for hydrogen and 


tungsten at ie = 573°K and Tp = 296°K. This value seems more in 


‘ 
keeping with the values measured in the present experiment for both 


nickel and gold with.nitrogen. 
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Figure 8. Variation of Rotational Accommodation Coefficient 
of Nickel, Silver and Gold with Surface Temperature. 
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The results of Schafer and Riggert (16) show that for gold and 
nitrogen Op reduces from 0.78 to 0.54 as b increases from 293°K to 
389°K. In the two methods they employ for determining the accommodation 
by the filament, the filament is heated by hot gas reflecting from a 
heated band in one case, and in the other it is ohmically heated. 

Thus, in the former case, energy is being given by the gas to the solid, 
while in the latter, the heated solid gives energy to a colder gas. 

Both Sasaki et al and Schafer and Riggert have used indirect methods 

to deduce Op « The accuracy of such methods is hard to assess. 

Results of both Knudsen (21) and Marsden (18)' show that the 
rotational and translational accommodation coefficients are equal. 
Knudsen found that the difference in both translational and thermal 
coefficients for rough and smooth surfaces was equal for monatomic and 
diatomic gas. Marsden obtained Op = 0.85 for nitrogen on silver. It 
should be noted that this result was obtained with the gas heated and 
the surface cold in contrast to the results obtained from room 
temperature gas and hot surfaces in the present experiment. 

The amount of theoretical work that has been done for predicting 
rotational accommodation coefficients is very limited. Feuer's (8) 
theory, based on a quantum mechanical model, states that oy should be 
less than Qe. This has been found to be generaily true in the present 


work. 


Oman's (7) theory, based on computer simulation of the collision 
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process, states that the energy exchange to the rotational mode is 
linearly proportional to the energy in the translational mode of the 
incident gas. The theory is based on a solid assumed to be a lattice 
of discreet atoms at 0°K. It has not been possible to check this 
theory because the translation energy of the incident gas was constant 
and the surface was at a higher temperature in the present experiments. 

An analysis of the expected accuracy of the present experiments is 
given below, following which some aspects of the method used for 
determining the rotational temperature of the gas are discussed. 
5.2 Estimation of Accuracy of Measured Rotational 

Accommodation Coefficient 


Rotational accommodation coefficient Op is given by equation 1.6: 


Ce on oe ime ieee ice §) 126 
R Re R. S R. 


Both Tp. and qe were measured by calibrated thermocouples. The error 
i 
is not expected to be more than +1°K. Experiments conducted with 


Us = T, = 295°K (room temperature) showed that the scatter of experimental 


R. 
; 


points from a least squares best fit line had a standard deviation of 
3.2% corresponding to extreme values of rotational temperature 
~ 


of T, +9°K. The slope of the least squares fit line gave the 


R R 
( iG 
values of Tp to be 295°K +5°K. The results of a typical run are 
r 
shown in the appendix. The calculated value of Tp was found to be 


r 
290°K for this data. 


At higher target temperatures, Tp » measured by subtracting two 
r 
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sets of line intensities, one of which is 1.5 times the other, shows 
considerably more scatter. The computer program given in the appendix 
prints out the rotational temperature after points lying more than a 
standard deviation away from the least squares fit line are removed. 

It also prints out the extreme rotational temperatures corresponding 

to the variation in slope due to standard deviation. The results gave 
a scatter in the rotational temperatures of less than +4%. 

Since the scatter and the line intensities are independent of 


temperature and the slope of the best fit line is proportional to 


R° 


V/Tp » the error in Tp » namely ATp » Will be proportional to T 
iG 4 iF if 


Hence, the scatter will cause ATp to be +0.04T, : 
r 


r 
The rotational accommodation coefficient can be written as 


t 
5 ee R. 
Gp = (Tp + ATp - Tp M/(T, - Tp ) = dp + 0.04 oy ne mae 
r r i j | ree 
S R. 
Set 


The limits of error given by this expression have been plotted on 
Figure 9 which shows the results for stainless steel. It is evident 
that the possible error increases beyond 10% for i below 440°K. ‘The 
error decreases as ie increases and Op decreases. 

It seems appropriate at this stage to examine the accuracy of the 
method used for determining the rotational temperature. The method 
was developed by Muntz (19). Its accuracy has been questioned by 


Ashkenas (22) and Harbour (23) among others. 
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Experiments by Marsden (18), Robben and Talbot (24) and Ashkenas 
(22) indicate that the measured rotational temperature is too high. 
The error increases with increase in density of the gas and decrease 
in its temperature. Ashkenas further suggests that the temperature 
depends on the number of peaks considered. 

Ashkenas concludes that above Tp = 300°K, the rotational temperature 
can be measured with less than 10°K error without applying corrections he 
has suggested fer density, temperature and the number of peaks. Based 


on a number of results, a statistically corrected equation of the form 


/T = A+B logo 5.2 


T Spectrometer true 
has been developed. The values of A and B are given in Figure 10 of 
his paper. At 300°K, for the number of peaks (10 to 13) and density 
‘ea ag gms per cubic cm.) used in the present work, A = 1.36 and 

B = 0.05. The equation 5.2 then becomes 

ih Rages 

1.002255 Das 


/T 


‘Spectrometer true 


it 


This shows that at 300°K, which is the lowest gas temperature used in 
the present experiments, the error is about 0.2%. For higher 
temperatures, the error is expected to be less than this. 

The discrepancy in rotational temperature has been attributed to 
secondary electrons produced, which do not necessarily follow the 


optical selection rule AK = +1. 
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Lillicrap (25) has proposed a theory in which the secondary 
electrons are taken into consideration by altering the equation put 


forth by Muntz (equation 17 of (19) ) to 


1(K')/2k'(G, + nG.) proportional to exp | - K'(K'+1)/T, 5.4 


where Gy is the same as [G] in Muntz' equation for primary electrons, 
G. is similar to Go for secondary electrons, and n is the number 
density of the gas. The above equation reduces to that of Muntz as 
n+0Q. Experimental data cited by Ashkenas (22) confirms the veracity 
of equation 5.4. 

Lillicrap has plotted a graph of line intensities normalized by 
dividing them by the highest line intensity, (1(3) at 80°K) against 
density. The plot shows that the use of his method gives good 
agreement for all densities. It also shows that equation given by 


y molecules per cubic centimeter. 


Muntz produces good results below 1 x 10 
The present experiments were conducted with a number density of 1.54 x 
nay molecules per cubic centimeter. 

Harbour (23) contends that there is a possibility of very high 
concentration of the secondary electrons. Experimental results of 
Dunn and Self (26) in parallel geometry collisionless flow are cited. 
It should be noted that their experimental set up was dissimilar to the 
one used during the present experiments. Mean energy of the secondaries 


is taken as 20 eV and the most probable energy as 0.6 eV. The secondaries 


are assumed to reach equilibrium temperature of 2 eV quickly. Harbour 
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argues that the electrostatic field of the ions reduces the radial 
velocity of the secondaries, and the field of the secondaries increases 
the radial velocity of the ions. This argument does not apply to the 
present experimental set up as the primary beam is not confined and 

there is no reason for ions to congregate near the beam centre. The 
field of the primary should stay negative. The extent to which Harbour's 
argument holds may be the reason the results improve at higher temperatures. 
In a geometry like the one presently used, the ions which inhibit the 
radial velocity of the secondaries should move out of the region of 
interest more quickly at higher temperatures because their own velocity 
would be very near that of the thermal agitation velocity of the gas 
molecules. 

Since the density in the present case is very much lower than those 
used in the experiments by most of the others mentioned, and the 
temperature equal to and above 300°K, the measured rotational temperatures 
are expected to be within the 0.047, scatter predicted by the analysis 


of experimental errors. 


5.3 Translational Accommodation Coefficient 
Figure 10 shows the variation of translational accommodation 
coefficient Oe with target surface temperature T.: for nickel, silver 
and gold. The translational accommodation coefficients are 0.5, 0.4 
and 0.2 at 400°K respectively for nickel, silver and gold. The variation 


of Oh with Me approximates to the equation 
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da. /dT = -13.7278 x 10° 
where wy is the ratio of the molecular mass of nitrogen to the atomic 
mass of the material of the surface. 

Among the works cited before, the results of Sasaki, Taku and 
Mitani (13) show that at T. = 791°K and 7, oto ae = 0.393 for 
nickel and nitrogen. Present results give Qt = 0.18 for the same 
conditions. The difference can probably be attributed to the difference 
in the condition of the test surfaces in the two experiments. Sasaki, 
Taku and Mitani calculated dy from the rotation of the torsion balance 
when a beam of molecules impinged on the test surface. The method used 
by Stickney (2), consisting of rotating of the balance shaft to keep 
the test surface in the same position, appears to be a better method 
of conducting this type of experiment. 

Sasaki and Mitani (14) have given the value of a, = 0.48 for 
hydrogen on tungsten at 573°K. Neither of the above mentioned works 
give the variation of Or with Ue 

Results of Schafer and Riggert (16) for gold and nitrogen show that 
Oy decreases from 0.89 to 0.86 as te increases from 293°K to 389°K. The 
values appear rather high. 

Marsden (18) has given de = 0.85 for silver and nitrogen, but his 
values are for a heated gas on a cold surface. 


The more interesting result of the present study is the behaviour 


of Oy with changes in surface temperature. There appears to be a 
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Figure 10. Variation of Translational Accommodation Coefficient 
of Nickel, Silver and Gold with Surface Temperature. 
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linear variation between ly and Tee the slope da,/dT being given by 
equations 5.5% 

Some data which support the behaviour of Oy predicted by equation 
9.5 iS given now. Carpenter, Humphries and Mair (15) have reported 
results which show the slope do/ dT. for nitrogen on tungsten to be 
Yan rime while equation 5.5 gives the slope as 3.38 x 10°" for the 
Same conditions. 

Results of Roach and Thomas (27), using the temperature jump- 
method, show the slope to be positive for neon and helium on tungsten, 
as predicted by the equation, though the magnitude of the slopes do not 
Show agreement. There is however, good agreement for argon, the slope 


4 


predicted by the equation for argon being - 0.47 x 10 © and the results 


giving a slope of - 0.38 x ln e 
Results of Thomas (28) using the low pressure method, show the 


: for mercury on platinum. The equation predicts 


the slope to be - 11.1 x ieee The results are for mercury which was 


Slope to be - 7.7 x 10 


not strenuously cleaned and Thomas states that the results are probably 

for mercury on an adsorbed film of mercury. The accommodation coefficient 
is probably high at higher temperatures due to evaporation of the 

adsorbed molecules. 

Equation 5.5 gives good agreement with results obtained in this 
experiment for nitrogen on a variety of metals, but when applied to other 
combinations such as helium on tungsten, it gives values of da,/dT, which 
are too large. It would appear that the constants in equation 5.5 are 


functions of parameters other than Megs Equation 5.6 is proposed to 
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replace equation 5.5. 
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This equation satisfies the results of the present experiments and 
agrees well with most of the experimental results cited previously. 
Since not many results are available in literature showing the variation 
of Oe with aes a better check on its validity could be made if it were 
to be integrated to give Oy directly as a function of ihe and other 

known parameters, such that 


Oy = (da,/dT. ) Te - Constant 57 


where the constant of integration is independent of Te The constant 
was determined by fitting equations 5.6 and 5.7 to experimental data 
given by Thomas (29), for helium, argon, and krypton on clean tungsten 


and bulk potassium. The resulting expression is given by equation 5.8 


a; 
Constant = = (61.19M-12525.2 YM + 28451.63)(32+Y2 (Y+1)(1.67-Y)M-M, ) 
M(184-Mc) 
- 5.6224M - 99.47 - 
° (62.07v - 5.6224M - 99.47 OM-10 85 /1+33.6 
(Ty, ~ 84.45viT + 4.7583M + 342.87) 5.8 


In this expression, the ratio of the specific heats were used 
to make it fit the results of the present work for diatomic nitrogen. 

The real test of empirical equations such as those given by 
equations 5.6 to 5.8 is how well they predict results for a wide range 
of gas-surface combinations and conditions, other than those used to 
determine the constant. This is discussed below for monatomic, diatomic, 


and polyatomic gases separately. 
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Monatomic gases: 
The table below shows some of the results predicted by the equations. 
These have been compared to experimental results given by Thomas(29). 


The difference between surface and gas temperatures have been taken 


to be 17°K. 


TABLE 3 


Predicted accommodation coefficients for monatomic gases 


Combination TT, °K a,(predicted) a(meas ured) 


‘ 
pa , cama 
i 


of 


The experimental results for helium on nickel are taken from Goodman's 
(4) paper. It appears from the table 3 that the equations predict 
good results for the gas-surface combinations given there. It was 
found, however, that for neon on tungsten they give Or varying from 
0.2 to 0.14 for a temperature range of 153° to 303°. These values 
are much higher than those given by Thomas (29) for the same conditions 
but are about half those predicted by Goodmans (4) original simple 
cubic lattice theory. 

It can be stated in general then, that the equations work 
well for a large range of gas-surface combinations for monatomic 
gases, neon excepted. 
Diatomic gases: 
As stated before the equations have been adjusted to give an excellent 
description of the results obtained here for nitrogen on gold, silver 
and nickel. Of the other gases tested with the equation, hydrogen 
on platinum gave values of Oe varying from 0.17 at 100° to 0.2 at 
600°K, the temperature of the surface being 17° higher. Though the 
trend of d- versus ie is opposite to that Slit by curves 7 and 9 of 
figure 2 of Thomas' (28) paper, the overall values compare well. Oxygen 
on platinum, on the other hand, gave values of Or which were lower than 
those given by Thomas (28), varying between 0.42 at 100°K to 0.22 at 
400° compared to about 0.5 at 400°. Weinberg and Merrill's (30) experiments 
have shown that platinum behaves like nickel when reflecting light gases 
such as helium and like silver when reflecting heavier gases such as 


argon and this may explain the discrepancy between measured and 
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predicted results in this case. It is unfortunate that dependable 


results for other metals are not available in literature. 


Polyatomic gases: 


The only gas tested with the equation in this category was carbon 
dioxide. The equation predicts Or on platinum to be 0.45 at 

300°K and 0.33 at 400°K, compared to 0.48 at 300°K and 0.36 at 400°K 
given by Thomas (28). At higher temperatures, however the experimental 


results are higher than those predicted. 


Obviously equations 5.6 to 5.8 need to be checked more 
thoroughly before their validity is established, but the results 
given here are deemed sufficient to support the proposition that 


Cy varies linearly with He according to equation 5.5. 
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9.4 Estimation of Accuracy of Measured Translational 
Accommodation Coefficient 


From theory, equation 4.17 gives 


| - | Mey, - 1) 4.17 

i 1 

The temperatures eg and Ty were measured by calibrated thermocouples 
] 

and the error in them is not expected to be more than +1°K. 


Experiments conducted with Ts = Ty = 295°K (room temperature) 
indicated that the scatter in (I, - 1,)2 was + 1%. Since there is 
ifr 


' 
to be any different the scatter in it 


LF 


no reason to expect (I, - I fe 

Z ] der 

r 

will be assumed to be + 1%. Scatter in (I, - I,) 
is calculated below, based on this. 

Letga 5 la )e eG 
. : =a 
yr = (I, I,) it and p= 7/ 


Then an error in p can be written as 


Ap = a a5 hr 5.9 
if 


where Aq ia = AY}, =+ 0.01 as noted before. 


Results of ms: typicalirun gave values of q = 163¢ units and r = 138 


units corresponding to Vr = 790°K. 


Substituting these values 


2 Pa 

163 163°x 138 _ 

Ap = —==x 0.01 + — x 0F0] =,0.03 5.10 
P1382 138¢x 138¢ 


T.3 net f | a gii\ a | 
esiquoocmiant batsidi fas yd bowessH stew apt hes T eowis: . 
AST net stom sd oF Batosqys Son ef my 10118 bn 
(swgsrsqnes moos) 3920S = {7 = oF rity bataubnoo atmambyegrd 
et siedd sonfe 8 + 26w ae ~ gf) ni sedtso2 of tedg beteotbat 
tt af s9ads22 ont Inatsttib one ad of atl ~ of) doaqxs of no2sst on : 
rity = QD ate ~ gi) mt sierihe’ tt 4 ecot Semuees ; 
‘ ett no. beead .woled bsxefuotesiah 
A gis tod +o 
NP om Ge gl) =" ; 
es neddiw sd mao q ni tore as nedT 
e2 1A < - R= gp 
-ovoted bston 26 10,0 += a\¥4 = o\ Pa srattw 
Ser = + brs ag tau ear = p to 2aulsy svsp muy [sotqya 6 Yo et fueeh ae 
MOR = Tot poibnoqzsrroa ettn - 
‘ lad seont pnitustiedue 
ate €0.0 = 10,0 x eta tal 


60 


Then hy Can be written as 


op = (p + Ap - WAT 7 - 1) = oy eA apre elie ye 
i j 
This variation is plotted on Figure 9. It is evident that the results 
fall within this probable error. 
It should be noted here that (I, - ty. is measured only once 
; 

for each set of experimental runs. Hence, the values of Oy at various 
lies will be affected to the same extent by the error in this variable in 
each run. 


For each run p can be written as p = q/r, where q is constant. 


Then iN) = a. Ar 9.12 
r 


Substituting values indicated before 


163°. 138 
Ap = Kea XU UT = 0 UO lico 
Pp 138° 1382 


Then Q- can be written as 


+ 0.0125/(T, Ay, 2) 5.13 


yp = Op 
Therefore values of Oh for each experimental run have a probable error 
+0.0125/(T. - 1), when their magnitudes at different T. are compared 
to each other. In other words, the slope da./dT. is more accurate 
than the absolute value of Oy. 

Throughout this work Oly and Op have been defined on the basis of 
temperatures. Such a definition assumes that both the incident and the 


reflected gas are in equilibrium. The justification for making such 
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an assumption is as follows. Equation 4.10 is derived from equation 
4.9 which depends on Maxwellian distribution of the molecules. It has 


been stated previously that a plot of In|I(K')/(K' + K'41)Ce](2 ) 
0 


against K'(K' + 1) will result in a straight line only if the gas has 

an equilibrium rotational energy distribution. Figure 7 shows such a 

plot and is typical of the results obtained during the experiments. It 
can be seen from the figure that a straight line is obtained, confirming 
equilibrium distribution of the rotational energy. This justifies the 

use of temperatures to define the rotational accommodation coefficient 

as in equation 1.6. With regard to the translational energy distribution, 
the incident gas molecules are in Maxwellian distribution. Hence, whether 
they are reflected diffusely or specularly the number reflected in 
various directions should still correspond to equilibrium distribution. 

The results of molecular beam experiments carried out by various authors 
indicate that the velocities of the reflected molecules also retain 
equilibrium distribution. For these reasons equation 1.5 defining 


Or based on temperatures is expected to be well justified. 


5.5 Surface Conditions 

The values of both translational and rotational accommodation 
coefficients are dependent on the surface conditions. Measurements 
of energy accommodation by molecular heat transfer methods have shown 
that adsorped gas layers can change the thermal accommodation by a 


factor ot iG. 
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Wachman (12) has reported 0.02 values of thermal accommodation 
coefficient of helium on clean tungsten, 0.185 for helium on tungsten 
with gas layers. No data is available indicating the effect of gas 
layers on the rotational coefficient. 

To obtain results for clean metals, efforts were made to conduct 
tests on surfaces which were (a) pure as far as material composition 
was concerned and (b) free from adsorped gas layers during the 
experiment. The composition of nickel, silver and gold sheets used 


for target material are given in Table 4, as supplied by the manufacturer. 


Table 4 
Typical Analyses of Surface Material Composition 


Serercola: Nickel 


Fine Silver 


Au 99.99% Ag 99.9% Ni 99.9% 

Copper .005% Cu .09% C .01% 

Silver .001% Pb .03% Mn .001% 
Fe .005% Fe .01% 


Ni 002% 
Zn 005% 


Cd .005% 


Al .005% 
Bi 01% 


Total Other 
.06% 
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The steps taken to keep the surface clean have been detailed in 
experimental procedure. All the experiments were conducted at high 
temperatures, and at these temperatures no gas layers are expected to 

be physically adsorbed. As far as chemical adsorption goes, Table 5, 
taken from Lewin (33), shows that at room temperature silver is inactive 


but nickel and gold are active. 


Table 5 


Chemical Adsorption 


where Y indicates chemical adsorption and N indicates none. 


Oil from the diffusion pump is the likely source of most of the 
hydrocarbons. Two cold traps on the roughing pump line as well as the 
length of the pump line are expected to prevent oi] vapours from getting 
into the system. Ho and 05 are formed by dissociation of water vapour 
in places where high temperatures are generated such as near the 
ionization gauge. Baking of the vacuum system should exclude most 
of the water vapour. On the basis of this discussion it is expected 


that the results are for clean metal surfaces. 
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CHAPTER VI 
CONCLUSIONS 


1) A technique for the simultaneous measurement of rotational and 
translational accommodation coefficients of nitrogen on a variety of 
solid surfaces has been developed. Reliable and repeatable results 


have been obtained with an accuracy of 


Oy Ls Or + 0.03/(T./Ty, a 1D) and Op = Op + 0.04(a, + Tele ~ TR) 


2) Values of dp were appreciably less than Chr » in agreement with 
theory, but the absolute values are much greater than might be expected 
from the few theoretical results that predict quantitative values. They 


also showed little variation with To 


3) Values of es on the other hand, varied with Te Results from 
three metals, nickel, silver and gold could be correlated by the equation 


a 


do/ dT. = i euicro x10) stec-004Znx deen 


Integrating this equation and adjusting the constant of integration 
using available experimental data, the resulting closed form solution 


for translational accommodation coefficient Cy» 
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not only gives an excellent description of the results obtained here, 
but also predicts correctly for a large range of gas-surface combinations 
and conditions. 

The proposition that the translational accommodation coefficient 


varies linearly with surface temperature appears to be valid. 


a) aks a, can be measured by this method for other gases and surfaces, 
it would be interesting to check the above equation for gas-surface 
combinations such as helium, and argon on various surfaces. One of 
the advantages of the present technique is that it places few 


restrictions on the material of the surface. 


5) High temperatures of the surface can be employed. Theoretical 
error analysis shows that the accuracy increases as the difference in 


temperature between the gas and the surface increases. 


6) Absolute pressure measurements are not required. This not only 
facilitates the performance of the experiments but also adds to the 


inherent accuracy. 
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APPENDIX 


The mechanism of excitation of the nitrogen molecules from the 
ground state Ny xls to an ionized state by the electron beam has been 
given in detail by Muntz (31). Only the highlights will be mentioned 
below. 

It has been stated in reference (19) that the limiting pressure 
at which a direct excitation to NoB“z State and subsequent emission can 
take place, unaffected by gas kinetic collisions, is about 470 x ee 
torr at room temperature. The pressure used in these experiments was 
5 orders of magnitude lower. 

In none of the experiments did the temperature of the gas exceed 
800°K. Hence, there was no significant excitation of the vibrational 
mode (19). Figure 11 shows schematically the excitation process at 
the top, the emission process in the middle, and the excitation emission 
path voltage levels at the bottom, and is reproduced from reference (19). 

The ionization, implying the carrying away of a secondary electron 
with a spin of 1/2, is the cause of change in multiplicity of 1 during 
the excitation from No xls to NoBoD. The spin coupling in the upper 
state is very weak, and it can be assumed that the doublet states are 
unresolved. The transition is assumed to belong to Hund's case (b) 
(32) and the applicable selection rule AJ = +], except for the spin 


change of 1/2 due to the removal of an electron. The rotational levels 


are designated K = J+1/2 and the selection rule, then AK = +1, gives 
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Figure 11. Mechanism of Fluorescence. 
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two paths of excitation, P and R, as shown in Figure 11 (top) corresponding 
to populating of energy level K' of the upper state from two levels of 
the ground state whose energy levels correspond to K' +e rvThus, as 


shown in Figure 11, the rotational energy level of N5BoE state, K' = 2 


I 
— 


is populated from the energy levels of the lower state No xls, Kits 
and 3 along R and P branches respectively. AK = 0 transitions are 
forbidden. 

The emission process is somewhat similar to the excitation process. 
The transition is from NoB“E to No x*5, The middle part of Figure 11 
Shows the transition schematically. There is no multiplicity change 
and the lower state also has a very weak spin coupling. The 5; to 
ey transition reduces the 5 to |S transition, if the doublets are 
unresolved. This again corresponds to Hund's case (b) (32) type, and 
the selection rule AK = +1, will lead to P and R branches of emission. 

A theoretical calculation of line intensities has been carried out 
by Muntz (19). This has not been reproduced here. One notable feature 
of the line intensities is that in the ground state No xls of the nitrogen 
molecule, the nuclear spin causes the odd rotational levels (K'' = 1,3.5, 
etc.) to have exactly half the population of the even levels (19). 
Selection rule AK = +1, populates the odd numbered levels in the upper 
state with twice as many molecules as the even numbered. This causes 


the peaks corresponding to K' = 1,3,5, etc. to be twice as intense as 


those corresponding to K' = 2,4,6, etc. 
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Only the R branch of emission is of interest as it is the only 
branch the spectrometer could resolve. The P branch forms a bandhead 
and turns back on itself, making it very difficult to resolve. 

Given below are the wave lengths at which the peaks of the R 
branch occur. These values correspond to rotational energy levels K' 
from 3 to 15. They have been taken from a paper by Child (20). The 
difference in wave lengths between successive peaks in Angstrom units 
is also given. This forms the criterion for working out the dimensions 


of the slits. 


Table 6 
Wavelengh of peaks in the R-branch 


Wavelength A Difference 
Fee ee oe Te 
07.533 
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As stated, the spectrometer could resolve only this band. By 
resolution, it is meant that when one of the peaks is centred over the 
exit slit, the neighbouring peaks do not influence its intensity. 

The resolving capacity of the set up was dependent mainly on 
three interrelated factors. The first was the inherent capacity of 
the spectrometer to resolve. This depends on the linear resolving 
power of the spectrometer and the widths of the entrance and exit slits. 
The details for calculating proper slit widths are given by Marsden (18) 
and these have been incorporated into the set up. The slit widths 
were a compromise between obtaining good resolution and intensity of 
light signal. Too large a width resulted in poor resolution and too 
narrow a width reduced the signal to below acceptable signal-noise 
ratio. The widths were set at 0.003 inches for entrance and 0.004 
inches for exit slits, as these gave the best results. 

The second factor affecting resolution was the speed at which the 
Spectrometer scanned. Too fast a scan resulted in poor resolution because 
of the time constants inherent in the coupled electrical measuring 
devices. Internal friction of the scanning mechanism prevented it 
from being run below a certain speed. 

The third factor is related to the second: the time constants of 
the measuring devices used, namely the electrometer and the chart recorder, 
and noise produced by the photomultiplier and the two above-mentioned 


instruments. The fluctuations in the output of the EMI 9502 S photo- 
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multiplier depends on the light intensity being measured. If, for 
example, the light’ intensity causes a nominal 1000 pulses per second, 
and this is averaged over periods of 1 second, then the root mean 
Square of the fluctuations is statistically equal to the reciprocal 
of the square root of 1000 or about +3%. The greater the line intensity, 
the higher the number of pulses produced and the lower will be the 
statistical fluctuation. To integrate out this noise, a 500 mFu 
Capacitance was inserted across the anode and the ground of the 
photomultiplier near the electrometer terminals. 

An integrating RC circuit was inserted between the electrometer 
and the chart recorder. It consisted of a 10 K potentiometer and a 
1 mF capacitance. The potentiometer resistance was adjusted until the 
noise from the electrometer was damped out. The electrometer was on 


10 x 10°18 


ampere scale. Its output was three volts. The chart 
recorder was adjusted to give full-scale deflection at two volts. It 
was offset to read about 1/10 of the full scale deflection when the 
electrometer was reading the dark current which had a value of 4 x eee 
amperes. Power supplied to run the spectrometer scanning mechanism 
motor was adjusted so that the scan between peaks a K' = 1 to 1/7 took 
about 10 minutes. 

A computer program was written to calculate the rotational tempera- 
ture. The data supplied were the line intensities at K' = 3 to 15, as 


given in Table 2. The following steps were carried out by the program 


which is attached. 
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Ik The temperature dependent factor [G] was calculated at Tp = 300°K 
for all values Of CKO) GCUni Lika / Gk? + Ki" + Disha were then 
calculated. (Zo a constant for each value of poe yee stored in the 
program. 0 

2. A least squares straight line was fitted to these points with 
K'(K' + 1) as the abscissa. From the slope of this line, Tp was 


calculated. If this T, differed from 300°K by more than 25°K, the 


R 
new temperature was used iteratively in step 1. 


3. The variance was calculated and all points lying outside this 

were thrown out. Step 1, 2 and 3 were repeated until all outliers 
were eliminated. Then the rotational temperature, extreme temperature 
corresponding to the standard deviation and the percentage deviation 
were printed. 

For purposes of checking the number of points used to define the 
line, the print out contained data points used, the points thrown out 
being printed as zeroes. A print out of the results for the case 
when the surface temperature is equal to room temperature is also 


attached. 
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DIMENSION A(30) ,B(40) ,¥(30), X(30), 
ST(10) ,@(30) , YCAL(30) 
READ(5.20)N,NUMBR 
FORMAT(2110) 
READ(5.22)(B(1) ,1=1,17) 
FORMAT(8F10.3) 
B(I) WAVELENGTHS OF R-BRANCH PEAKS 
DO 2 I=1,17 

B( I)=3908.2/B(1) 
PELL TOPS: 

B(1)=B(1)*B(T) 
X(I)=1*(1I+1) 
CONTINUE 
B(I)=NEW/NEWNOT RAISED TO 4 
X(1)=K'(K'+1) 
KOUNT=0 
READ(5.21)(A(1) ,1=3,N) 
FORMAT(8F10.3) 
A(1) MEASURED LINE INTENSITIES STARTING AT K'=3 
READ(5.112) TEMP 

FORMAT (F10. 3) 
WRITE(6,111)TEMP 
FORMAT(//T8, ‘SURFACE TEMPERATURE, DEGREES KELVIN=',F10.3) 
VCR) =02 
G(1)=0 
L=0 
M=0 
T(1)=300. 
B(19)=2.00065*1.4384 
B(19)=BV1' '*HC/K' 
BV1I=ROTATIONAL CONSTANT IN GROUND STATE 
H=PLANCK'S CONSTANT 
C=SPEED OF LIGHT 
K=BOLTZMANN'S CONSTANT 


K=0 

DO 4 I=3,N 
IF(A(1) .£Q.0.)GOTO70 
IF(X(I).EQ.0.)GOT060 

SEL EIN AEY P(-2.*B(19)*I+1)/T(1)) 
G(1)=G(1)+I*EXP(2.*B(19)*1I/T(1)) 
G(1)=G(1)7 (22 *1r1) 

EUVSHORTOTEC THER 

Y(1)=ALOG10(A(1I))-G(1)-ALOG10(2.*1) 

GOTO60 

¥(1)=0 

M=M+1 
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WRITE(6,45)G(1I) ,¥(I) ,X(1) ,A(1) 
FORMAT(/ ,4F10.3) 

G(I)=(G) OF MUNTZ'S EQUATION 
Y(I)=LOG 1(K')/(K'+K' '+1)(G)B(1) 
CONTINUE 

SUMY=0. 

SUMX=0. 

SUMXY=0. 

SUMX X=0. 

DO 5 I=3,N 

IF(A(I) .EQ.0.)GOTO5 
IF(X(1I).EQ.0.)GOTO5 
SUMY=SUMY+Y(T) 

SUMX=SUMX+X(T) 
SUMXY=SUMXY+(X(1)*Y(T)) 

SUMX X=SUMX X+(X(1)*X(1)) 

CONTINUE 
WRITE(6.49 )SUMY ,SUMX , SUMXY , SUMXX 
FORMAT (//,4F14.3) 

SUMY X=SUMY*SUMX/ (N-M-L-2) 
SUMXP=SUMX* SUMX/(N-M-L-2) 
SLOPE=( SUMXY-SUMY X) / (SUMX X-SUMXP ) 
CONST=SUMY-(SUMX* SLOPE) 
CONST=CONST/(N-M-L-2) 
WRITE(6.51)SLOPE ,CONST 
FORMAT(2F10.6) 


T(2)=2.00065*1 .4384/(2.302585*SLOPE) 


TC2)=ABS (12) ) 

TB) SARS G2) =O) 

Po =aheZ) 

K=K+l 

IF(K.GI.7)GOTO11 
WRITE(6.40)K,T(1) 
FORMAT(//T8,110,'APPOROX, TEMP.= 
M=0 

IF(T(3) .GT.25.)GOTO3 


SUM=0. 

DO 6 I=3,N 

YCAL(1)=SLOPE* 1*( I+1) +CONST 
IF(A(1).EQ. a )GOTO6 

rF(X( ) .EQ.0. ) GOTO6 

SUM=SUM+(Y (I ie YCALGL))**2 

CONTINUE 


SIGMA=2.*(SUM/(N-M-L-2)**.5 
SIGMA=2.*SIGMA 
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SIGMA=VARIANCE 
WRITE(6,61)SUM, SIGMA 
FORMAT(2F 10.3) 

WRITE(6,75) (YCAL(I) ,1=3,15) 
FORMAT( F10.3) 

J=N 

C=SUMY 

DO 9 I=3,N 
IF(X(I).EQ.0.)GOTO9 

IF(A(I) .EQ.0.) GOTO9 
IF(ABS(Y(I)-YCAL(I)) .1.T.SIGMA) GOTO9 


C=SUMY-Y(I) 
SUMX=SUMX-X(1) 
SUMX Y=SUMXY - X( 
SUMXX=SUMXX =X ( 
NL )=0% 
X(1)=0s 
L=L+] 
CONTINUE 
IF(C.EQ.SUMY)GOTO11 
SUMY=C 
GOTO 10 
WRITE(6.25)T(1) 
FORMAT(//T8,'ROTATIONAL TEMPERATURE=' ,F10.3) 
XBAR=SUMX/(J-2) 
XSQ=0 
DO) ta el=3,N 
IF(A(I).EQ.0.)GOTO12 
IF(X(1).EQ.0.)GOTO12 
XSQ=XSQ+( X(1I)-XBAR)**2 
CONTINUE 
SIGMA=SUM/(J-2) 
SIGSP=SIGMA/XSQ 
STDEV=SQRT(ABS(SIGSP) ) 
STDEV=STANDARD DEVIATION 
T(4)=ABS((T(1)*SLOPE) /(SLOPE+STDEV) ) 
T(5)=ABS((T(1)*SLOPE) /(SLOPE-STDEV) ) 

) 


Yeh) 
1)*X(I) 


WRITE(6,30)T(4),T(5) 

FORMAT(//T8,'EXTREME TEMP.=',2F10.3 
T(4)=(T(4) -T(1)) /T(1)*100. 
POS)=CTS) HT CO) 7 11) 2100; 

WRITE(6,31)T(4) ,1(5) 

FORMAT(//T8,'PERCENTAGE TEMP. DEVIATION=' ,2F10.3) 
KOUNT=KOUNT+1 

DO 00 s1=15 17 

X(1)=1*(1I+1) 

IF( KOUNT.LE.NUMBER )GOT052 
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